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Abstract 

We describe electrical transport in ideal single-layer graphene at zero applied bias. There is 
a crossover from collisionless transport at frequencies larger than ksT/h (T is the temperature) 
to collision-dominated transport at lower frequencies. The d.c. conductivity is computed by 
the solution of a quantum Boltzmann equation. Due to a logarithmic singularity in the collinear 
scattering amplitude (a consequence of relativistic dispersion in two dimensions) quasi-particles and 
-holes moving in the same direction tend to an effective equilibrium distribution whose parameters 
depend on the direction of motion. This property allows us to find the non-equilibrium distribution 
functions and the quantum critical conductivity exactly to leading order in 1/| log(a)| where a is 
the coupling constant characterizing the Coulomb interactions. 
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I. INTRODUCTION 



Despite the intense experimental and theoretical interest in the electronic properties of 
graphene Q], there has been relatively little progress in measuring and understanding the 
role of electron-electron interactions. However, the recent ability to grow ultrahigh mobility, 
suspended, single layer graphene [2I, 0] promises that the situation may well change in the 
near future. 

This paper will examine the role of electron-electron interactions in an infinite sample of 
single layer graphene without impurities. We will also restrict our attention to the undoped 
case, so that the chemical potential is at the node of the massless Dirac spectrum. Our 
results can be extended to include a non-zero chemical potential and a dilute concentration 
of impurities: this was discussed recently in Ref. Q for a low-frequency 'hydro dynamic' 
regime, and additional results will appear in forthcoming work. 

The key to understanding electron-electron interactions in clean, undoped graphene is 
the fact that it is a nearly 'quantum critical' system with marginally irrelevant Coulomb 
interactions This implies that the inelastic, electron-electron scattering rate 

is of order ksT/h, where T is the absolute temperature, and there is a crossover from 
'hydro dynamic' to 'coUisionless' transport as the measurement frequency {u) is increased 



past the scattering rate 
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These two regimes are captured in the following limiting 



forms for the frequency dependence of the electrical conductivity, a, 
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1 



ln(a(T))| 



where a(T) is a temperature-dependent, dimensionless 'fine-structure constant', which con- 
trols the strength of the electron-electron interactions (defined more precisely in Section HTl). 
and A is a cutoff energy scale of the order of the electronic bandwidth. The high frequency 
result above (the 'coUisionless' regime) was obtained in Refs. lo], 0. The leading term is the 
conductivity of 4 species of free m;.sles= Di.ac Herbut al Q abo obtained the 

coefficient of the subleading [ln(A//ki;)]~^ term. The low frequency result in Eq. (11. ip . which 
is the collision-dominated hydrodynamic regime, is the primary new result of this paper. At 
asymptotically low temperatures we have (see Eq. (12. 4p ) 



a{T) 



ln(A/T)' 



:i.2) 



the resulting logarithmic increase of a with decreasing T is similar to those of quantum 
critical systems in their upper-critical dimension ISj, and the inelastic scattering rate of the 
carriers is of order [ksT /K)a^{T). 
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Related results have been obtained recently by Kashuba IJ] in a preprint which appeared 
while our paper was being completed. 

We also note that our results are obtained in the context of a solution of the quantum 
Boltzmann equation. Going beyond the Boltzmann approximation, and in a system with 
perfect momentum conservation, we have to consider potentially singular hydrodynamic 
"long-time tails" in a mode-coupling theory 15|], which could modify the low frequency 
behavior of the conductivity. Such effects are however innocuous here, because the long- 
range Coulomb interaction suppresses density fluctuations. 



II. RENORMALIZATION GROUP ANALYSIS 

Here, and in the remainder of the paper, we set = /c^ = 1. 

We begin introducing the low-energy theory for graphene, and reviewing its renormal- 
ization group (RG) properties. The theory is expressed in terms of = 4 species of 
two-component Dirac fermions (a = 1 . . . A^) and the Euclidean partition function 



V^^VAr exp(-5) 



a=l 



— + 2eA^(x, r) + iv^pa'' — + i-A^ 
or \oxc 
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The functional integral is over fields defined in two spatial dimensions x = (x, y) and imagi- 
nary time r, cr^'^ are Pauli matrices acting on the sublattice space of the honeycomb lattice, 
and v^p is the bare Fermi velocity. The scalar potential, A^, mediates the e^/ (e|x|) Coulomb 
interaction between the electrons, where e = {ea + £b) /2 is the dielectric constant for a 
graphene sheet confined between two dielectrica with dielectric constants Ea and Eb, respec- 
tively. We have also introduced an non-fiuctuating external vector potential A = {Ax,Ay) 
as a source field: this allows us to extract the electrical current. 

The renormalization group properties of Z have been discussed elsewhere aaaa 
The fermion field undergoes a wavefunction renormalization, the charge e remains 
unrenormalized, and the velocity vp renormalizes to larger values with decreasing energy 
scale. For the velocity renormalization, we have the RG equation 

^ = fic^W, (2.1) 
where the running "fine-structure constant" is 



a = , (2.2) 

EVp 
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and the function /(a) = a/A in the perturbative regime of small a. We can re-express these 
results in terms of the RG equation for the dimensionless coupling a 

Notice that a scales to small values at small energies, and this is what facilitates the transport 
analysis of this paper. It has been shown that a = is the only fixedpoint in an analysis 
which, in the large N limit, also remains valid for large values of a j^, H]- 

We are only interested here in observables related to the electrical current, and so we will 
not need the explicit form of the wavefunction renormalization. The current is obtained by 
taking a functional derivative with respect to A, and this is protected by gauge invariance 
to have the same form when expressed in terms of either the bare or renormalized quanti- 
ties [l^, which we will use explicitly in fl3.11|3.12l) below. For two dimensional graphene 
this implies that the scaling dimension of the conductivity is exactly zero and is unaffected 
by wavefunction renormalizations. This result can also be obtained explicitly by exploring 
charge conservation of the system along with the related Ward identity |9| and holds to 
arbitrary order in perturbation theory. 

We are interested here in the collision-dominated transport regime, where the character- 
istic energy of excitations is ksT. We thus use the RG equation to scale down from some 
high energy cutoff scale. A, to a scale ksT. Integrating Eq. (12.31) over this regime, we obtain 

^ ' l + (a74)ln(A/T) ~ ln(A/T) ' ^'^ 

where a" is the bare value dependent upon Vp. Son 15] has also examined the structure 
of the RG flow at strong coupling in the large N limit; he finds that there is a significant 
intermediate energy scale over which 



a{T) ~ . (2.5) 

Both Eqs. (12. 4p and (12.51) predict a slow flow with decreasing temperature towards weak 
coupling. We can also use a(T) to obtain a T-dependent velocity 



VFiT) 



a(T) 



Yln(A/r) 



(2.6) 



We also note that the leading order flow in a in Eq. (12. 3p represents an exchange- 
correlation effect. Ordinary screening effects are formally higher order, and can be accounted 
for in the random-phase approximation by the mapping 



16 



171 



4 



III. COLLISION-DOMINATED TRANSPORT 



After initially renormalizing down to a scale T as described in Eq. ([ITj), we can now inves- 
tigate the transport quantities in the renormalized theory. So all subsequent appearances of 
the field the velocity vp-, and the coupling a implicitly refer to the T-dependent renor- 
malized quantities obtained as described in Section [Tll We will not explicitly write-out this 
T dependence. 

Our formulation of the transport properties of the renormalized theory of weakly- 
interacting massless Dirac fermions closely follows that presented in Ref. [hI. This previous 
work considered massless Dirac fermions interacting with a weak statistical interaction due 
to a Chern-Simons term, and here we only need to replace the Chern-Simons term by a 
Coulomb interaction. The transport analysis is easiest in the real-time operator formulation 
with the Hamiltonian 

H = H^ + H^ 

Ho = I [vp-^i {-la'd,) -^a] (3.1) 



(3.2) 



with the Coulomb interaction 



^(q) = (3.3) 

I T.I 

and a = 1, . . . ,N labeling the "flavors" of fermions (A^ = 4 in graphene, accounting for 2 
valleys and 2 spin projections). Even though we compute our results specifically for the 
Coulomb interactions (13. 3p . the formalism carries through in exactly the same manner for 
arbitrary isotropic two body potentials. 

The simplest formulation of the transport equations is in a basis which diagonalizes the 
Hamiltonian Hq. To do this, we first express \1/ in its Fourier components 

^.(x,*)- /T^^'f'!! (3.4) 



(27r)2 \^C2.(k,() 

and then perform a unitary transformation from the Fourier mode operators (ci^, C2a) to 
(7+a,7-a): 

Cla{k) = -^(7+a(k)+7-a(k)) 

cUk) = -|^(7+a(k)-7-a(k)). (3.5) 

We have introduced here a notational convention that we shall find quite useful in the 
following: as k is a two-dimensional momentum, we can define the complex number K by 

K = kx + iky where k = {k^, ky) (3.6) 
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and k 
result 



|k| = \K\. Expressing the Hamiltonian Hq in terms of 7±, we obtain the simple 

d^k 



A, a 



;XvFk-fl^{k)-fXaO<-), 



(3.7) 



where the sum over A extends over +, — . 

Let us also express the interaction Hamiltonian Hi in terms of the jxa- 

(Pki <fk2 <Pq 



Hi 



A1A2A3A4 



(27r)2 (27r)2 (27r) 



X T: 



Ai A2A3A4 \ 



(ki, ka, q)7Lb(ki + q)7l3a(k2 - q)7A2a(k2)7Aib(ki) (3.8) 



where 



1 A2A3A4 (ki,k2,q) 



V(q) 



1 + A1A4 



{Kl + Q*)Ki 
|ki + q|/ci 



1 + A2A3- 



2 - q|fc2 



(3.9) 



Finally, we also express the electrical current, obtained by taking a functional derivative 
of the action with respect to A, in terms of the 7±. For the case of a spatially independent 
current (which is the only case of interest here), the result can be written as 



J = J/ + J 



with 



and 



Xa 



(fk Ak 
(2^T 



7L(k)7Aa(k) , 



-levp 



(Pk (z X k) 



7t(k)7-„(k)-7la(k)7+a(k)' 



(3.10) 



(3.11) 



(3.12) 



(27r)2 k 

where z, a unit vector orthogonal to the x, y plane. J/ measures the current carried by 
motion of the quasiparticles and quasiholes — notice the A prefactor, indicating that these 
excitations have opposite charges. The operator J// creates a quasiparticle-quasihole pair, 
and in the low frequency limit of interest here we may neglect J//. Similarly to the problems 
studied in Refs. [lO, [lH , a current carying state with holes and electrons moving in opposite 
directions is consistent with a vanishing total momentum. Thus a finite conductivity does 
not require the total momentum of the problem to relax. This is the physical reason why at 
the particle hole symmetric point, i.e., at vanishing deviation of the chemical potential from 
the Dirac point, the d.c. conductivity is finite even in the absence of momentum relaxing 
impurities. However, as we will see below at finite deviation from particle hole symmetry a 
driving electric field always excites the system into a state with finite momentum. The 
latter cannot decay which entails an infinite d.c. conductivity, in accordance with the 
hydrodynamic analysis 
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We can now write down the collisionless transport equations for the excitations. As a 
first step, we define the distribution functions 

/A(k,t) = (7L(k,t)7Aa(k,t)\. (3.13) 



where there is no sum over a on the rhs, and we assume the distribution functions to be the 
same for all valleys and spins. In equilibrium, i.e., in the absence of external perturbations, 
these are related to the Fermi function 

/4k, t) = r(^i^fc) = -^^^,3^L___ 

/.(k,t) = n-v,k) = (3.14) 

where we temporarily allow for a finite chemical potential /i. 

Then to first order, in the presence of an external electric field E, we find the simple 
equations 

| + eE-|^)/A(k,t)=0. (3.15) 



It is a simple matter to solve fl3.15p in linear response. First we parameterize the change in 
fx from its equilibrium value by [isj] 

/A(k,u;) = 2Ti5{uj)f\\vFk) + e^^^^f\XvFk){l - f\XvFk))gx{k,u), (3.16) 

where we have performed a Fourier transform in time to frequencies, cu, and introduced 
the unknown function gx{k,uj). At the particle hole symmetric point (/i = 0), an applied 
electric field generates an deviations in the distribution functions having opposite sign for 
quasiparticles and quasiholes, 

gxik,uj) = Xg{k,uj). (3.17) 

This reflects the fact that there is an increased number of quasiholes and quasiparticles 
moving parallel and antiparallel to field, respectively. As quasiparticles and -holes have 
opposite charges, their electrical currents are equal, while their net momenta have opposite 
signs. 

Inserting (13.161) into (I3.15p . we obtain a simple solution for the function g 

9ik,co)= ""^'l (3.18) 
[-IUJ + rj) 

where t] is a. positive infinitesimal. Inserting this result into fl3.13p and (13.111) . we obtain the 
conductivity 

^ ' E{uj) {-tuj + 7]) J {2ny k^ \ dk J 

NkeT \n2 

h {-ihw + T]) ^ ' 
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where, in the last equation, we have re-inserted factors of h and fc^. Note that all factors 
of vp have cancelled out: this is a consequence of the conductivity having scaling dimension 
(i — 2 (where d is the spatial dimensionality), and being independent of the dynamic critical 
exponent z. So in this free electron approximation, the real part of the low frequency a is 
a delta function at u; = with weight of order ksT. 

Including interband transitions the real part of the conductivity becomes 



e2 



Rea(uj) = —N-n 
h 



[/i;BTln2](5(M + ^tanh (-^ 



(3.20) 



with high frequency limit Recx (tu ^ ksT /K) c^Nti / {8h). In the collisionless regime this 
constant value remains the leading contribution to the conductivity even if one includes the 
electron-electron Coulomb interaction j^, [l^]. Next order corrections are of the form 

Nti 

Re a (cu > kBT/h) = -— [1 + O {a {u))] , (3.21) 

ll o 

where a{uj) ^ 4/ ln{A/hw) is the renormalized frequency- dependent 'fine-structure constant' 
for ksT <^ hLj <^ A. For = 4, this yields the result given in the upper row of Eq. (11. II) . 
Thus, in the collisionless regime, interactions only lead to very small changes of the conduc- 
tivity. In the following subsections, we will discuss the opposite, collision dominated hmit, 
hw -C ksT, and determine how collisions broaden the delta function of Eq. fl3.20p to a Drude 
peak. 



A. Quantum Boltzmann equation 



We now include collision terms on the right hand side of fl3.15p . We can determine these 
terms by application of Fermi's golden rule ll|], or by the explicit derivation presented in 
the Appendix: 



| + .E.A)A(k.«) 



(27r) f d'^ki d\ 



(27r)2 (27r)2 

6{k - A;i - |k + q| + |ki - q|)i?i(k, ki, q){/A(k, t)/_A(ki, t) [1 - ^(k + q, t)] 

X [1 - /_,(ki - q, t)] - [1 - /,(k, t)] [1 - /_,(ki, t)]A(k + q, t)/_,(ki - q, t) } 
6{k + A;i - |k + q| - |ki - q|)i?2(k, ki, q) |/A(k, t)/,(ki, t) [1 - /^(k + q, t)] (3.22) 

X [1 - /,(ki - q, t)] - [1 - A(k, t)] [1 - ^ik^, t)]A(k + q, t)Mk^ - q, t) } | . 
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where 

i?i(k,ki,q) = 4(|T+__+(k,ki,q)-T+_+_(k,ki,-k-q+ki)f 

+ (iV - l)|T+__+(k,ki,q)|' + (iV - l)|r+_+_(k,ki, -k - q + ki)|') , 

/I, ,2 

i?2(k,ki,q) = 4 ( -|T++++(k,ki,q) - T++++(k, ki, ki - k - q) | 

+ (iV-l)|T++++(k,ki,q)f) , (3.23) 



which are illustrated in Fig. [H The terms proportional to -Ri represent collisions between op- 
positely charged particles, while those proportional to R2 are collisions between like charges. 
There are also processes where a particle-hole pair is created: as in Refs. 10, lll|, these 



can be dropped because they have vanishing phase space upon imposition of the energy 
conservation constraint with dispersion Sk = vpk 




+,i +,i +,i +,i +,j +,j 



FIG. 1: Illustration of the Golden rule diagrams entering the collision term. The diagrams (a) 
describe scattering of oppositely charged particles corresponding to the term while the diagrams 
(b) describe scattering of like particles corresponding to the term R2. Note that the vertex preserves 
the flavor a = but not the particle/hole nature A = it. The factor 1/2 of the first diagram 
accounts for the symmetry factor associated with having two indistinguishable particles in the final 
state. 



We now proceed to the linearization of (13.231) by inserting the parametrization (I3.16P and 
find 
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-iuogx{k,uo)-XvF/T) k_ (27r) f (fki (fq 



(e^Ffc/T + i)(e-^ffc/T + i)^ j (27r)2(27r) 

5(A; - /ci - |k + q| + |ki - q|)i?i(k, ki, q, 

'"k N ki ,, , (k + q) . , (ki — q) . 

-gx{k,uj) + —g^x{ki,uj} - — — -^A(|k + q|, a;) - r^-A(|ki - q|, 

k ki |k + q| |ki — q| 

Sjk + fci - |k + q| - |ki - q|)i?2(k, ki, q) 

^ (^e-^pk/T i)(g-t>Ffci/T _|_ x)(efF|k+q|/T + i)(e^F|ki-q|/r _j_ x) 

k /, N ki (k + q) I , (ki - q) 

-gx{k,uj) + — 5(A(/i;i,u;) - |]^_^ 5'A(|k + q|,^) - jj^^ _ ^A(|ki - q|,t^) 



UJ 



(3.24) 



The remainder of this paper is focused on the solution of the hnearized transport equation 
in Eq. (13.241) for the function g. It is useful at this point to recall some crucial mathematical 
properties of such transport equations, reviewed, e.g., by Ziman 19|] and Arnold et al. 18|. 
We can view the right hand side of Eq. (I3.24p as a linear operator, the so-called collision 
operator C, acting on the function (k//c)(y'(fc); we drop the implicit uj dependence because C 
is independent of tu. A key property of C is that it is Hermitian with respect to the natural 
inner product 

{gM =Y, f 4i^9iAk)g2Ak). (3.25) 



(2vr) 



This Hermiticity follows 



from symmetry properties of Ri and R2 under exchanges be- 
tween incoming and outgoing momenta, which are very similar to those used in establishing 
Boltzmann's H-theorem. 

Related to the above properties of the collision operator, we can introduce a functional 
Q[g], such that Eq. (13.241) is equivalent to finding its stationary point 



ma] 
Sg 



0. 



(3.26) 



Specializing to the particle- hole symmetric case, cf. Eq. (13.171) . the explicit form of the 
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functional is 



_ (2^) f ^^'^ ^^'^i 
~ 8^ y (27r)2 (27r)2 (27r)2 



r 



X 



X 



^(fc - fci - |k + q| + |ki - q\)Ri{k, ki, q) 

^k , ki , (k + q) , , (ki - q) , 

-g{k, uj) - —g{k,, uj) - ^(^(Ik + q|, ^) + |i^^_q| ^(|ki " q| 

5{k + fci - |k + q| - |ki - q|)i?2(k, ki, q) 

■^-vpk/T I'^i^f.-vpki/T ]_)(^gi,_F.|k+q|/T _^ (^gt;^ |ki -q|/T _^ 

'k , ki . fk 



-g[k,u)) + -^^(fi^i, 



■^7(|k + q|, 



(ki - q) 



+ 



k + q| ' ^" ^ |ki — q| 

dPk g{k,uj)[—iujg{k,uj)/2 — VF/T] 
(27r)2 (e''ffc/T + i)(e-^Ffc/T + i) • 



^(|ki - q| 



(3.27) 



B. Translational invariance and momentum conservation 



The translational invariance of the system immediately implies the presence of a zero 
mode of the operator C, which corresponds to the shift of the distribution functions arising 
from changing to a linearly moving reference frame. The corresponding deviation gx{k,uj) 
has the form 

gx{k, uj) = C{uj)iJtrik) = C{uj)k, (3.28) 

which is easily seen to annihilate the left hand side of fl3.24p due to momentum conservation. 
Note that this zero mode of the Boltzmann operator is orthogonal to any modes of the form 
f l3.17p which are the only ones that can be excited by an electric field at particle hole 
symmetry. This again expresses the fact that current and momentum are independent of 
each other at this special point. However, away from the Dirac point, or if a thermal gradient 
is applied instead of an electric field, the zero mode tptr will be excited by the driving field 
which leads to a diverging d.c. response in clean systems. This will be discussed in more 
detail in a forthcoming publication. 

For the following we restrict to electrical conductivity at the particle-hole symmetric 
point where the above zero mode is irrelevant. 



C. Collinear limit 



In the previous analysis of a quantum Boltzmann equation for massless Dirac fermions 
in two dimensions it was noted that the phase space for scattering of particles was 
logarithmically divergent in the collinear limit. For the interaction considered in that paper, 
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the coUinear scattering cross-section vanished, and so this singular phase space density had 
no important consequences. The coUinear scattering does not vanish for the present Coulomb 
interaction, and so we need to consider this logarithmic divergence seriously. 

The physical origin of the divergent coUinear scattering is related to the linear dispersion 
which implies that quasiparticles or -holes moving in the same direction share the same 
group velocity, independent of their energies. This leads to a diverging duration of collisions 
of nearly coUinear particles, which is enhanced due to the low space dimensionality. To the 
extent that coUinear scattering is very strong, and if we consider small enough frequencies, 
we may expect that quasiparticles and holes that move in the same direction in the plane 
will establish a pseudo-equilibrium characterized by an effective chemical potential and an 
effective temperature which will however depend on the direction of motion. 

In linear response the deviations of these effective parameters from the equilibrium values 
II and T have to vary with k//c • E for symmetry reasons. Further, the effective temperature 
shift is easily shown to be identical to the mode iptr discussed above, and thus, it is thus 
ruled out at the particle hole symmetric point /i — 0. The remaining dominant mode of the 
function g will correspond to an effective shift in chemical potential which translates into 

gxik,u;) = Ciu;)ij,{k) ^ ^AC(a;), (3.29) 

where the prefactor has been chosen so as to make C{u;) dimensionless. With this Ansatz, 
which will be confirmed below, it simply remains to determine the prefactor C (uj) , yielding 
the leading term in the non-equilibrium distribution. Note that the effective chemical po- 
tential shift ranges between ±C hvpeE/T depending on the direction. Comparing this to 
the temperature allows us to estimate the threshold electric field strength, eE^^ — T'^/hvp, 
below which non-linear effects should remain small. 

Let us now review in more detail, how the above physical picture arises in the formalism of 
the Boltzmann equation. The occurrence of a logarithmic divergence can be seen by allowing 
the incoming and outgoing momenta to be nearly coUinear. Without loss of generality, we 
choose k = {k,0) with A; > 0. Also, we write ki = {ki,k±), q = {q,q±) with k± and q± 
small. The divergence in the phase space density of the coUision term proportional to R2 
occurs when ki > 0, k + q > and ki — q>0. Likewise, for scattering of oppositely charged 
particles a divergence occurs when their k vectors are anticoUinear which ensures coUinear 
group velocities since Vk = Ai'^k/A;. In this regime, the argument of the energy conservation 
delta function of the particle-particle scattering term can be written as 

kl ql (fcx - q±r 



^ + ^1 — |k -I- q| — |ki — q| 



2ki 2{k + q) 2{ki - q) 
{k + ki) 



2{k + q){ki-q) 



{q±-Cik±){q±-C2k±), (3.30) 



where Ci,2 depend upon k, ki, and q, and are the roots of a quadratic equation which 
are defined by the expressions above. Then, the phase space density for the R2 term is 
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proportional to 

dk± 4{k + q){ki - q) 
^(^ + A:i)|Ci-C2| 



j dk±dq±S{k + fci — |k + q| — |ki — q|) = J 



ki{k + q){ki — q) f dk 



\ki 



(3.31) 



The logarithmic divergence as fc^ — > is now evident. This divergence is clearly a conse- 
quence of the linear dispersion of the fermions, and the above analysis also makes it clear 



that it is special to two dimensions. As discussed in Ref. Il8| for a similar divergence in a 
different problem, we expect that this divergence is cutoff by higher-order self energy cor- 
rections to the fermions. Such self-energy corrections appear at order a in the perturbation 
theory, and so the important range of the k± integral is between T/vp and Ta/vp. So we 
may approximate 18j 

dh , 

' i21n(l/a), (3.32) 



\k±. 

and set k± = q± = elsewhere to obtain the leading contribution to the collision integral 
in the limit a ^ 0. Proceeding in this manner, the part of C on the right-hand side of 
Eq. (I3.24P proportional to R2, which we denote C2, becomes 



2n'^VF k Jq J _^ g2 y 

R2 [gjk, uj) + gjki, uj) - g{k + g, ^) - gjki - g, uj)] 

^ i)(g-i>Ffci/T + i)(e^F(fc+g)/T + i)(e''F(fci-g)/T _|_ 1) ' ^ ' > 

Consonant with our discussion earlier in this subsection, a key property of the above ex- 
pression for C2 was noted by Kashuba 1^: the function g = constant is an eigenvector 
of C2 with zero eigenvalue. The same is also easily seen to apply to the portion Ci of C 
which is proportional to Ri. Indeed, this is just the direction-specific chemical potential 
shift in Eq. (13.291) . which naturally is a zero mode for collinear scattering, since it maintains 
a pseudo-equilibrium among particles moving in the same direction. 

Going beyond the collinear limit, we conclude that there is an eigenvalue of C which is 
not proportional to ln(l/a) in the limit of small a; the corresponding eigenvector is given 
by a constant g{k) up to corrections of order [ln(l/a)]~^. 

The solution of the Boltzmann equation in Eq. (I3.24p requires that we obtain the operator 
C~^, and the results above allow us to constrain its form in the limit ln(l/a) ^ 1. Let 
be the eigenvectors of C with eigenvalues A^. Then 

C-i ^ ^ MM , (3.34) 



and in the limit of large ln(l/a), C^^ is dominated l4] by the eigenvector whose eigenvalue 
is not proportional to ln(l/a). Note that it is quite remarkable that in this limit we can 
solve the Boltzmann equation essentially exactly. 
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D. Results 



From the reasoning in the previous subsection, we conclude that up to corrections of 
order [ln(l/a)]~^, we can choose g to be of the form 



(3.35) 



We insert this parameterization into the functional Q[g] in Eq. (13.271) : the solution of the 
stationarity condition in Eq. (I3.26P is then equivalent to requiring the vanishing of the 
derivative with respect to C. We numerically evaluated the integrals in Eq. (13.271) using an 
elliptic co-ordinate system to solve the energy conservation constraint HI], and obtained 



Q[9] 



1 ln2 

T 17 



t l'^) duo] 



with K = 3. 646 for the physical case iV = 4. From the stationarity condition we then obtain 

1 



C{u;) 



(3.36) 



—i{uj/T) + Ka"^ 

The conductivity can be obtained from C{u!) by combining Eqs. (13.111) . (I3.16P and (13.351) : 

-2 A^A;Brin2 



a{uj) 



(3.37) 



h —ihu + KkBTa^ ' 

where we have re-inserted factors of h and ks- Notice that the conductivity depends only 
upon a(T), while all other factors of vp{T) cancel. Notice also the connection to the free 
particle result in Eq. (I3.19P — the only difference is that the infinitesimal rj has been replaced 
by the inelastic relaxation energy nksTa"^. 



IV. CONCLUSIONS 

We conclude by briefly noting the conditions under which our main results for conduc- 
tivity in Eqs. (I3.37p . (12. 4p . and (12. 5p may be observed in transport measurements. The key 
requirement is that ksTa^ be the largest infrared energy scale which quenches the ideal 
Dirac fermion behavior. Thus, the sample size should be larger than the inelastic scattering 
length £ee ~ hvp/ {kBTa^). Similarly, the elastic mean-free path from impurity scattering 
should be larger than 4e, too. Finally, particle- hole symmetry is also required, and so the 
bias voltage should be smaller than fc^Ta^. 

It is possible to extend our analysis to include all the additional perturbations noted in 



the previous paragraph, with a treatment of disorder effects following that of Ref. l20|. When 
these perturbations are weak (compared to ksTa"^), then in the collision-dominated regime, 
a general hydrodynamic analysis is possible: this was presented recently in Ref. Q. Also, 
in this regime the analysis of the Boltzmann equation greatly simplifies if the interactions 
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are weak enough to ensure a strong logarithmic divergence in the colhnear channeL The 
latter establishes pseudoequilibrium along different directions if the inelastic scattering time 
remains the shortest relevant time scale in the problem. Otherwise, a full analysis of the 
modified quantum Boltzmann equation is required. These aspects will be discussed in future 
work. 
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APPENDIX A: DERIVATION OF THE QUANTUM KINETIC EQUATION 

An alternative derivation of the quantum kinetic equation can be carried out in the 
framework of closed time contour ordered perturbation theory, as explicited in Chapter 9 of 
Ref. [2l|. The problem we consider falls into the generic class of a system describing particles 
interacting via a distance-dependent density- density interaction. Thus the Hamiltonian is 
of the form 

H = Ho + ^j d\d\'V{v, r')p(r)p(r') , (Al) 

where p(r) denotes the particle density at spatial point r. The following considerations 
are completely generic and apply to any system which falls into the class of Hamiltonians 

21 



presented in Eq. flAip . The starting point of our discussion is given by Eq.(9-7a) of Ref. 
which has to be generalized to incorporate a possible matrix structure of the Green's function 
(in our case the Green's function lives in spinor space within a structure due to the N spin 
and valley species, and thus has a 2N x 2N structure) 

[9t- VRf/(R,T)Vk]G<(k,c^;R,T) = - G<(k, cu; R, T)S>(k, cu; R, T) 

+ G>(k,cu;R,T)S<(k,cu;R,T) (A2) 

where 

S>^<(r,t;R,r) ^ -i^ I dRrff V(R + r/2-R-r/2) V(R-r/2-R + f/2) x 



X G<f(-r 



'75 
- G>'< 



r, -t; R, T) [G>^(r, t; R, T)G>^<(f , t; R, T) 
R + f/2 - R + r/2, t- R, T)G>^<(R + r/2 - R + r/2, t; R,(A^) 
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in the Bom approximation (note: double indices are summed over). The corresponding 
self-energy diagrams are the RPA-type contribution and the maximally crossed diagram, 
see also Ref. 2l|. In a next step, following the treatment in Kadanoff and Baym, we find 
that the Fourier transform with respect to the relative coordinates (which corresponds to 
the mixed Wigner transform) of Eq. flA3|) reads (note that in the following we drop the 
dependence on the centre of mass coordinate R) 



SXr(k,w;T) = , " \ .0 , ^o (27r)^,5(k + ki - ka - kg) x 



X 6{uJ + UJi — UJ2 — iOs) X 

X [Vik - k^)V{k - k2)G<f (ki, cui)G>^(k2, a;2)G>^<(k3, u,) 

- V(k-k2)V(k-k3)G<f(k,c.i)G>^(k2,u;2)G>^<(k3,a;3)] • (A4) 

Until now all the formulae are completely generic and not specific to graphene. In order 
to make connection to the problem of graphene we note that the Green's function of the 
spinors ^ is related to the Green's function of the 7 through 

G<'>(k,cu) = f/kV'>(k,cu)f/k 

where the unitary matrix t/j^^ according to Eq. (13. 5p is given by 

Furthermore we note, that the summation over spin and valley indices only affects the RPA- 
like diagram, which thus receives a prefactor and the resulting matrix equation is an 
equation, whose indices only carry over the 2x2 spinor. This allows to rewrite Eq. (]A4ll as 

,>^<, , f d^ki duoi d^'k2 duj2 d^k^ duo^ 3 




X 5{uJ + UOi — 002 — ^3) X 

X [NV{k-k2)V{k-k2) X 

X (f/,><'>(ki,^0f/kj^, (f^k>>'<(k2,^2)f/k.)„, (f^k>>'^(k3,^3)f/k3),, 

- y(k-k2)F(k-k3) 

X (f/k-'><'>(ki,^i)f/kj^, (f/k->>'<(k2,^2)f/k.)„^ (f/k->>'<(k3,^3)f/k30^) 

Accounting for the fact that the operators 7 describe sharp quasiparticles the lesser and 
greater Green's functions are given by 

^< ,(k, u- T) = 2ti5{lo - eA(k, T))h{k, T)5x,x' 

and 

^> ,(k, a;; T) = 2'n5{uj - tx{k, T))[l- A(k, T)] 5a,a' , 
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where we assumed the distribution function of the quasiparticles to have no off-diagonal 
components, which is justified to hnear order in the potential gradient. We can formulate 
the kinetic equation for the diagonal part of the distribution function as 

[dr - VR^7(R,T)Vk] ^(k,T) = - /^(k,T) (f/kS>(k,cu = e,{\^)-T)U^^) 

+ [1 - /^(k,T)] (t/k5]<(k,a; = 6^(k); T)^/^-^) JA6) 

or equivalently 

[dr - VRt/(R, T) Vk] ^(k, T) = - /^(k, T)a>^(k, u = e^(k); T) 

+ [l-^(k,r)]a<^(k,a; = 6^(k);r). (A7) 

Exploiting the form of the lesser and greater Green's functions we can rewrite the self- 
energies as (note that /x, in contrast to the other double indices, is not summed over here 
and subsequently) 

<(k,c. = e,(k);T) = y ^ — ^ — ^ — (2.) 5(k + k, - k. - ks) 

X (5(e^(k) +U1-LJ2- CJ3) X 

X [NV{k - k2)V(k - k2){2nf6{uJi - ex{ki))6{L02 - exA^2))6{LO, - e^M) 

X MA,A(k3,ki)MAA,(ki,k3)Af,a,(k,k2)MA,^(k2,k)/A(ki,T)(l-/A,(k2,T))(l-/A,(k3,r)) 

- V{k - k2)V{k - k3)(27r)35(u;i - ex{ki))d{u2 - e,., (k2))5(cu3 - ^xM) 

X MAA,(ki,k3)MA,A(k2,ki)M^Ai(k,k2)rA,^(k3,k)/;,(ki,T)(l - /A,(k2,T))(l - /A,(k3,r))] 

(A8) 

and 

< /, n\ -D rr\ f ^^ki dLUi d:^k2 duj2 d'^ks du^ 3 

.<,(k,a; = e,(k);R,r) = y ^ — ^ — ^ — (2.) 5(k + k, - ka - k3) 

X 5(e^(k) + u;i - c^2 - ^^3) x 

X [NV{k - k2)\/(ki - k2){2'Kf6{uJi - eA(ki))5(o^2 - ex,{k2))S{uj^ - ex,{k^)) 

X M;,,;,(k3, ki)M,;,, (kj, k3)M^,, (k, k2)M,,^(k2, k) (1 - /,(ki, T))/,, (ks, T)A,(k3, T) 

- V{k - k2)V{k - k^){2'Kf5{uji - ex{k^))5{uj2 - exA'^^m^^ - ex,{k^)) 

X MAA,(ki, k3)MA,A(k2, ki)M^Ai(k, k2)M;,,^(k3, k)(l - ^(ki, r))/Ai(k2, r)/A.(k3, T)] , 

(A9) 

where the shorthand notation 

1 / 1 _L 1 _ K*Ki \ 1 / K*K \ 
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was introduced. The connection with the matrix elements defined in Eq. (13.91) can be easily 
established and reads 

T,A,A.A3(k,ki,q) = ^V(-q)A/r,,3(k + q,k)MA,A.(ki - q,ki) . (AlO) 

Using Eq. flA7p and preforming a sequence of transformations finally yields 

[aT-VR?7(R,T)Vk]^(k,T) = — / \ ^/ (/ifc + Afci - Ai|k + q| - A^lk^ - q|) 

X [Ary(-q)y(-q)MA,A(ki -q,ki)MAA2(ki,ki -q)M^Ai(k,q + k)MA,^(q + k,k) 

- V(-q)V(k - ki + q)MAA,(ki, ki - q)MA,A(q + k, ki)M^Ai (k, q + k)MA,^(ki - q, k)] 

[(1 - ^(k,r))(i - /,(ki,T))/,,(q + k,r)A,(ki - q,r) 

- ^(k,T)/,(ki,T)(l-A,(q + k,r))(l-A,(ki-q,r))] . (All) 

In a next step we will make connection to the Golden rule result of the main text. It is 
straightforward to see that using Eq. flAlOl) we can rewrite the above expression to yield 

[dr - VrC/(R, T)Vk] Uik, T) = — —J^—^5{f^k + \h - Ai|k + q| - A2|ki - q|) 

X 4 [iV|T^,,,Ai(k, ki, q)p - T^AA.Ai(k, ki, q)T;,;,^,^(k, ki, ki - k - q)] 
[(1 - Uk, r))(l - /A(ki, T))/a, (q + k, r)/A, (ki - q, T) 

- /,(k,T)/A(ki,T)(l - /A,(q + k,r))(l - /A,(ki - q,r))] (A12) 

Energy and momentum conservation restricts the valid combinations of particles and holes 
scattering, see Ref . [ill . , which simplifies the above expression. Applying all these sim- 
plifications and shifting the variables appropriately we obtain 

6{k - A;i - |k + q| + |ki - q|):Ri|^(k,t)/_^(ki,t)[l - ^(k + q,t)][l - /_,,(ki - q,t)] 

- [1 - f,{k, t)] [1 - /_^(ki, t)]^(k + q, t)/„^(ki - q, t) } 
6{k + A;i - |k + q| - |ki - q\)R,[f^{k, t)f^{k,, t) [1 - /^(k + q, t)] [1 - /^(ki - q, t)] 

- [1 - Uk, t)] [1 - Uk^,t)]Uk + q, t)Uk^ - q, t) } | , (A13) 

where 

R, = 4Ar(|T+__+(k,ki,q)p + |r+_+_(k,ki,ki-k-q)p) 
-4T+__+(k, ki, q)n + (k, ki, ki - k - q) 

-4r;_+_(k, ki, ki - k - q)T+__+(k, ki, q) (AM) 



18 



and 

= 4iV|T++++(k,ki,q)|2 

-4T++++(k, ki, q)n+++(k, ki, ki - k - q) . (A15) 

Performing the appropriate shifts allows to write 

= 4(iV-l)|T+__+(k,ki,q)|2 + 4(iV-l)|T+_+_(k,ki,ki-k-q)|2 

+4|T+__+(k, ki, q) - T+_+_(k, ki, ki - k - q)|2 (A16) 

and 

R2 = 4(iV-l)|T++++(k,ki,q)|2 

+2|T++++(k, ki, q) - T++++(k, ki, ki - k - q)|2 , (A17) 

which establishes the equivalence of Fermi's Golden rule and the Keldysh treatment, see 
Eq. I^Ml- 
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